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Abstract. We consider parabolic systems with nonlinear dynamic boundary 
conditions, for which we give a rigorous derivation. Then, we give them sev- 
eral physical interpretations which includes an interpretation for the porous- 
medium equation, and for certain reaction-diffusion systems that occur in 
mathematical biology and ecology. We devise several strategies which imply 
(uniform) L p and L°° estimates on the solutions for the initial value problems 
considered. 

1. Introduction 

In this article, we consider the following system of quasilinear parabolic equations 

(1.1) d tUl - A (Ai ( Ui )) + fi (x, t, Tt) = 0, in tt X (0, oo) , 

for i = 1, ...,m, where u — (u\, u m ) , tt is a bounded domain in M. N , N > 1, 
with sufficiently smooth boundary V := dtt (which is at least of class C 2 ), for some 
given functions Ai and fi. Denote by N m = {1, ...,m} and consider two mutually 
disjoint (possibly empty) subsets I, mi J m C N m such that /,„ U J m = N m . Equation 
(|l.lj) is subject to the following set of boundary conditions 

(1.2) d n u t + hi (x, t, it) — 0, on T x (0, oo) , i e I m 
and 

(1.3) Sidtm + d n {Ai (ui)) + g.i (x,t, it) = 0, on T x (0, oo) , i e J m , 

for some given functions gi and hi. Here Si > for i £ J m , and we may assume, 
without loss of generality, that Si = 0, for i G I m . The boundary conditions in 
(|1.2j) - (|1.3j) may be also mixed, that is, the boundary T may consists of two disjoint 
open subsets Ti and T2 on which the boundary conditions may be either of Dirichlet 
type or of the form (|1.2I) and (|1.3p . respectively. Finally, the model (|l.ip - (|1.3l) could 
be also generalized be letting the reaction terms depend on advection, by allowing 
the diffusion rates depend also on x and t, or in other various ways. As usual, we 
equip the system (ll.l[) - (|1.3[) with the initial conditions 

(1.4) Ui\ t=Q = u i0 in tt, Ui\ t=Q = v i0 on T, i e N m , 

where in general, we may have Mio|r 7^ v iOi * £ N m (i.e., if u^o is well-defined in the 
trace sense). 
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We aim to give some results which allow to deduce L°°-estimates for solutions 
of (|l.lj) - (|1.4p assuming that some sort of energy estimate is apriori known in L p - 
norm for some finite p. The main tool will be an iterative argument following a 
well-known Alikakos-Moscr technique combined with a suitable form of Gronwall's 
inequality, which can then be used to prove that the L p -L°° smoothing property 
holds for any solutions of the non- degenerate parabolic system (|1.1[) - (|1.4[) (e.g., at 
least when ai (u^ :— A i (ui) satisfies (jl.5[) below). In order to deal with the full 
degenerate case (|1.1[) - (I1.4D (at least in the case when a* (ui) = \ui\ Pz , pi > 0), we 
employ DeGiorgi's truncation method to prove the L p -L°° smoothing property. The 
precise statements of these results can be found in Section 3, see Theorems 13.11 and 
13.21 A rigorous derivation and physical interpretation of the system (|1.1[) - (|1.4[) shall 
be given below in Section 2. 

Why is it important to establish apriori (possibly, uniform in time) L°°-estimates 
from some given L p -estimate? To better give an idea of our larger scope let us take 
a look at some history for problems of the form (ll.ll) - (|1.4p . Problems such as 
(|1.4[) have already been investigated in a number of papers [16j H3 El EI] assuming 
that diffusion rates ai (ui) :— A i (ui) satisfy 

(1.5) ai (ui) > di > 0, for all m € R and i 6 N m . 

For instance, Constantin and Escher (161 117] show that unique (classical) maximal 
solutions exist in some Bessel potential spaces under suitable assumptions on the 
nonlinearities fi, gi and hi. Such results also enable the authors to investigate other 
qualitative properties concerning global existence and blow-up phenomena (see, also 
Ej). These results are also improved by Meyries [31], still in the non-degenerate 
case (|1.5[) . by assuming more general boundary conditions (by also incorporating 
surface diffusion in (|1.3p ), and by requiring that the functions fi ( u ) , hi (it) , gi (it) 
are dissipative in a certain sense. However, none of these contributions deal with 
the degenerate case for equation (jl.ip . that is, when ai (ui) is allowed to have a 
polynomial degeneracy at zero for some (if not all) i £ N m ; for instance, one can 
take 

(1.6) at (ui) = \u % \ Pl , pi > 0. 

Moreover, it is well-known in the scalar case m = 1, that when at least one of 
the source terms, the bulk nonlinear term f\ or the boundary term g\ is present 
in (|l.ip - (|1.2p . conditions can be derived on their growth rates which imply either 
the global existence of solutions or blow-up in finite time [3T]. Namely in the non- 
degenerate case, for A, /i 6 {0, ±1} with max {A, p} = 1, fi (s) := — A|s| ri_1 s and 
gi (s) := — p \s\ T2 1 s, solutions of 

(1.7) d t u - vAu + fi (u) = hi (x) , in x (0, +oo), 
subject to the dynamic condition 

(1.8) d t u + ubd a u + gi (u) = h 2 (x) , on T x (0, oo) , 

are globally well-defined, for every given (sufficiently smooth) initial data (|1.4p . if 
r\Ti > 1 and Ari + pr2 > 0. Furthermore, |21j shows that if we further restrict 
the growths of r x ,r 2 so that r x < (N + 2) / (N - 2) and r 2 < N/ (N - 2), then the 
global solutions are also bounded. These restrictions can be eventually removed 
and more general conditions on /i, <?i can be deduced (see, e.g., |23j). On the other 
hand, if A = 0, fj, = 1, then some solutions blowup in finite time with blowup 
occurring in the _L°°-norm at a rate (t — T*) ' r2 ^ , for some additional conditions 
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on uq and r 2 . In the same way, when /j, — and A = 1 , then some solutions blowup 
in finite time with a blowup rate which depends on n and uo (see [3]). In the case 
when both \i G K and A > are nonzero, blowup may still occur for superlinear 
growth of fi and any growth of g (see [23j). The occurrence of blow up phenomena 
is closely related to the blowup problem for the ordinary differential equation 

u t + H (u) = 0, 

where cither H = j\ or H = g\ . More precisely, it is easy to see that solutions of the 
ODE are spatially homogeneous solutions of cither equation (| 1 . T[) or (|1.8[) . and so 
if these solutions blowup in finite time so do the solutions of (|1.7p - (|1.8[) . It is worth 
mentioning that in [8] a criterion for the global existence of a (classical) maximal 
solution (on some interval [0,t+)) to (11.7|) - (|1.8j) is established using a variation of 
parameter formula. In particular, it is shown that if t+ < oo then necessarily we 
must have 

limsup||u(t)|| ioo(n) = oo. 

t — 

Therefore, it appears that in order to deduce global existence of classical solutions 
to systems of the form (|1.1[) - (|1.4[) . (|1.5[) . it is generally required that we should 
deduce bounds on the solutions in L°°-norm (see [3j] also). 

Finally, the L v -L°° smoothing property also becomes an essential tool in attractor 
theory where it can be used to establish the existence of an absorbing set in L°°- 
norm if this property can be deduced easily in L p -norm for some finite p (in many 
applications in physics and mechanics, p is equal to either 1 or 2). Recall that a 
subset BcH, where % is a topological space endowed with a given metric, is called 
absorbing if the orbits corresponding to bounded sets V of initial data enter into B 
after a certain time (which may depend on the set V) and will stay there forever. 
Moreover, we note that in order to study the long term behavior of the parabolic 
system (|1.1[) - (|1.2[) . if the absorbing property holds in L°°-norm, the growth rate of 
the nonlinearities fi, gi and hi with respect to Ui becomes nonessential for further 
investigations of attractors. Indeed, the absorbing property can be also established 
in higher-order VF S:P -norms with relative case provided that it is known in L°°-norm. 
For the application of this property to attractor theory for parabolic equations of 
the form (|1.7[> . (|1.8I) . see [22] [53], where explicit dimension estimates for the global 
attractor for (|1.7[) - (|1.8[) are obtained. 

The main goal of this paper is to deduce sufficiently general conditions on the 
diffusions and sources in (|l.l[) - (jl.3[l . which would prevent blowup of any solution 
in the L°°-norm, and show that the parabolic system under consideration is dissi- 
pative in a suitable sense. We outline the plan of the paper, as follows. In Section 
2, we give the full derivation of systems of the form (|l.l[) - (jl.3p . and give physical 
interpretations to the dynamic boundary condition (|1.2[) for the porous-medium 
equation, and some models in ecology. In Section 3, after we introduce some nota- 
tions and preliminary facts, we give the statements of our main results and some 
further applications. Finally, in Section 4 we provide the full proofs of these results. 

2. Derivation and interpretation 

Let r C R^ -1 consists of two disjoint open subsets Ti and T2, each T^r^ is a 
S'-null subset of T and r = Ti U T 2 with Ti C T. We shall only give the derivation 
in the case of the scalar equation 

(2.1) d t u - div (a (u) Vti) + / (u) = h Y (x) , 
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equipped with (nonlinear) dynamic boundary conditions 

(2.2) d t u + a (u) Vu • n + g (u) = /12 (x) , 
on Fi, and Dirichlet boundary conditions 

(2.3) u lF2 = 0. 

We can easily extend our arguments to systems as well (see below). Equations 
(|2.1[) - (I2.3[) are also subject to the initial conditions 

(2.4) Uu—o = uq in i7, Uu—o — v<j on T. 

Standard derivations of the porous medium equation always use the principle 
"amount of fluid in equals amount of fluid out" over a region f2, occupied by either 
a liquid or gas, and is based on the fact that this fluid diffuses from locations of 
higher to those of lower pressure. In the traditional approach, the porous medium 
equation is assumed to hold in the region f2 and then the boundary conditions 
are appended later. There are three standard boundary conditions that specify 
the density on the boundary of f2; they are Dirichlet and Neumann- Robin type 
of boundary conditions (see, e.g., [35J ). Dynamic boundary conditions for porous 
medium equations seem to have appeared before in different contexts [5J EH [Ml [35] • 
For instance, [6] deals with the modelling of the rain water infiltration through the 
soil above an aquifer in regimes where there is runoff at the ground surface. In 
general, all rain water infiltrates into the soil, but if the rainfall event is particularly 
intense, the maximum draining capacity of the soil is exceeded. In this case, dynamic 
boundary conditions (see (|2.8p below) are needed to describe the saturation of layers 
near the ground surface (cf. also [HI [36]). Porous-medium like systems (|2.I[) - (|2.3I) 
can also be found as part of (larger) coupled systems of partial differential equations 
(such as, (|l.l[l - (|1.4[) ) describing the vertical movement of water and salt in a domain 
splitted in two parts: a water reservoir and a saturated porous medium below it, 
in which a continuous extraction of fresh water takes place (for instance, by the 
roots of mangroves) [24) . Such problems are formulated in terms of equations for 
the salt concentration and the water flow in the porous medium, with a dynamic 
boundary condition which connects both subdomains. Finally, dynamic boundary 
conditions similar to (|2.2[) also appear in certain classes of parabolic equations 
with boundary hysteresis (see, e.g., [351 Section 4] and the references therein). For 
some applications of dynamic boundary conditions for physiologically structured 
populations with diffusion we refer the reader, for instance, to [20] . 

For all the phenomena of the kind discussed here, the method of introducing dy- 
namic boundary conditions seems ad hoc. It would be more natural if such bound- 
ary conditions could be derived in the context of energy balance and constitutive 
laws. Moreover, the usual derivation of the porous medium equation with standard 
boundary conditions does not show how to model, for instance, a water source, 
which is located on the boundary of tt. To this end, we shall rethink the usual 
derivation of the porous medium equation (12.11) . by making essential connections 
between the differential equation (|2.ip and the boundary conditions (|2.2[) - (I2.3I) . and, 
thus, try to convince the reader that our new perspective is more natural than the 
traditional way. Let p(x,t) denote the pressure of fluid at x £ Q and time t > 0. 
Consider the mass of fluid in an element of volume V given by 

a (x) u (x, t) dx, 
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where a (x) > defines the porosity of medium at the point x € f2. Similarly, we 
define 

J /3(x)v(x,t) dS 

as the mass of fluid across the surface T, where (3 (x) is such that 

r 3 := {x E ri : (x) > 0} 

is a set of positive measure and CTi. In what follows, we shall take T2 = for 
the sake of exposition, so that Ti = T. The flux JJ (x, t), at which the fluid moves 
across a surface element S with normal n, is given by 

/ I(x,t) ■ ndS. 
Js 

Suppose now there is a source on the boundary T to be represented by a function 

* = x, u, Vw). 

The amount of fluid leaving the region is still given by J r JJ (a;, t) ■ ndS, but the 
amount of fluid leaving into the region must also take into account the action of 
the source VP on T. It is worth pointing out that, in practice, when rainfall only 
partially infiltrates the soil, the water will accumulate on the ground surface T% as 
the surface layer becomes saturated; hence, necessarily, $ ^ on T. We use the 
measure space (fi, dp) which we redefine as (f2, dx)®(T, dS). Then the conservation 
of fluid in f2 takes the form 

(2.5) d t (j a (x) pdx + J (x) V dS^j + J I ■ ndS 

= j Edx+ j ^>dS, 
Jn Jr 

where S = S {x, t, u) denotes any volume source density function. Notice that 
equation (|2.5I) must also account for a term like 

J P(x) vdS, 

due to the presence of the source density ^ at T. Assuming that J is sufficiently 
smooth and applying the divergence theorem in (|2.5p . we deduce 

(2.6) dt( I a (x) udx + I /3 (x) udS) + [ div (J) dx 

\Jq Jr J Jn 

= j Edx+ [ ^dS. 
Jn Jr 

Assuming that the density functions u, v are also differentiable with respect to 
t > and since (|2 .6[) holds for any subdomain f2 Q the usual argument yields 
the following differential equation 

(2.7) d t (a (x) u (x, t)) = -div (J (1, t)) + 5 (x, t, u (x, t)) , x 6 fi, t > 0. 
Then, from (|2.6|) the boundary condition becomes 

J [d t (/3 (x) u(x,t))--$ (x, t, u (x, t) , Vw (x, t))} dS = 0, t > 0, 

which clearly holds if 

(2.8) d t (/3 (x) u (x, t)) - * (x, t, u (x, t) , Vu (x, t)) = 0, for x E T, t > 0. 
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Darcy's law states that the flux J depends on the pressure gradient so it takes the 
form 

(2.9) J (x, t) = (x, t) Vp (x,t), x en, t> 0, 

v 

where v > is the viscosity of the fluid and K. (x) defines the permeability of the 
porous medium. Finally, if one also makes the assumption that the pressure p is 
described by an equation of state involving the density, p = b (u), then substituting 
the appropriate quantities in (12.7[) . we obtain the following porous medium equation 
(2.10) 

dt (a (x) u (x, t)) = — div (/C (x) a (u (x, 1)) Vu (x, 1)) + S (x, t, u (x, t)) , in f2, t > 0, 
v 

where we have set a (t) = tb (t) . The function b that relates the density to pressure 
is, in general, monotone and, in fact, strictly increasing in many applications of the 
type of fluid being considered in the literature. 

Finally, let us now focus on the boundary condition (|2.8p . We now show that a 
quite large class of boundary conditions for equation (|2.1j) can be written in this 
way for various choices of VP. We emphasize that in this formulation the boundary 
conditions arise naturally in the formulation of the problem. Suppose first 1^ = Ti 
(i.e., P(x) > a.e. in I\) and j3 G C 1 (Fi) . Choosing \fr = (i.e., no source is 
located at Ti), so that dtu = on Ti = T, therefore 

(2.11) u (x, t) = u {x) , 

for x £ Ti and t > 0, where uq is the initial condition associated with equation 
(|2.ip . Thus, we obtain a Dirichlet boundary condition for u. In order to derive an 
inhomogeneous Neumann boundary condition, we suppose that ^ only depends on t. 
Then, if u is sufficiently regular, we have dtu(x, t) — (1//3 (x)) ^ (t) on Ti, for any t > 
0. Hence, if Ti is smooth enough as well, we have dt (Vw) = V (dtu) — 7 (x) "J (t) 
on Q x (0,+oo), for some 7(2:) e R N . This entails that X7u(x,t) = H(x,t) holds 
for (a;, t) G Ti x (0, +00) and some smooth function H e Mr. Therefore, we have 

(2.12) Vit (x, t) ■ n = H (ar, t) ■ n, (x,t) G Ti x (0, +00). 

To obtain a Robin boundary condition, we set * (x, t, u, Vu) = /3 (x) e Cr q (t) , for 
some C G R, where r is defined as the parameter describing the line I which passes 
through x and contains n such that r > at all points on I n £1 which are close to 
x. We thus have dtu (x, t) = e Cr q (t) which implies 

V (d t u (x, t)) ■ n (x, t) = d t (Vu (x, t) ■ n) (2, t) = Ce Cr q (t) , 

for (x,t) G Ti x (0, +00). Therefore, we infer 

9t (Vu (x, • n) (x, t) — C9 t u (x, t) = 0, on I\ x (0, +00), 

so that 

(2.13) Vu (x, t) ■ n - Cu(x, t) = j (x) , (x, f) G Ti x (0, +00). 

We have thus recovered the most common boundary conditions. On the other hand, 
in order to model a water source placed on the boundary Ti, we may assume that 
W only depends nonlinearly on the flux Vu • n across the boundary as well as on a 
nonlinear source g (u) , which may represent effects of reaction or absorption at Ti. 
That is, we let 

\& (x, u, Vu) = —a (u) Vu • n — g (u) — h 2 (x) . 
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Then, the resulting boundary condition for equation (|2.1I) becomes 

(2.14) (3 (x) d t u (x, t) + a(u (x, t)) Vu (x,t)-n + g (u (x, t)) = h 2 (x) , 

for (x, t) £ T\ x (0, +oo). This is a reasonably general condition which contains the 
usual (homogeneous) ones along with the so-called dynamic boundary condition 
when T3 = Ti = T. When ponding or surface runoff occurs at the surface Ti, 
it also includes the dynamic boundary condition contained in Filo-Luckhaus [19] . 
Finally, we notice that /3 (x) may be such that T3 ^ I\, so that we can also cover 
the case where the boundary conditions (|2.14j) arc dynamic only on a part of the 
boundary of T\. If the source '5 also depends on v in the tangential coordinates of 
the boundary, i.e., "J = ^ (x, u, Vu, Vr^ w) , then we can also model local diffusion 
in (|2.14[) by incorporating the elliptic Laplace-Beltrami operator Ar x (or other 
nonlinear differential operators) on the manifold Ti C M.^ -1 . 

We will now give a physical interpretation of the effect of a water source on the 
patch T±, at least in some cases. A similar approach was also used in the derivation 
of heat and wave equations (see, e.g., [25]). We will mainly focus on the following 
boundary condition 

(2.15) d t u + a (u) Vu ■ n = 0, on T t x (0, 00) . 

We work in an infinitesimal region on the boundary. Choose a point x £ T± and 
let B K (x) be a ball of radius n > about x. Since T\ is regular, we can choose 
a coordinate system for B K (x) (~l O so that the boundary of B K (x) n O in the 
transformed coordinate system is flat, x is mapped tox — (xi, X2, xn-i, 0), that 
is, the boundary T, at least locally near x lies on the hyperplane xn — 0. Then the 
outward unit normal n to T at x is the unit vector in the direction of e at which we 
will denote by r. Then, locally near x, (|2.15l) becomes 

(2.16) dtu + d r (A{u)) = 0, (r,t) E (0,r ) x (0,t ), 

for some sufficiently small positive constants ro, to. Observe that (12.16|) resembles 
nothing more than a scalar conservation law, where we have set A (u) = J Q a (t) dt. 
Equation (|2.16|) . subject to initial condition u (r, 0) = Vq (r) , r G (0, r ) , possesses 
interesting types of fundamental solutions such as, travelling waves that describe 
the movement of a mass of fluid in the direction of the unit normal n G M. N , and 
source-type solutions starting from a finite mass concentrated at a single point of 
space, say, vq (r) = CSq (r) , C > 0. 

In the latter case, explicit self-similar solutions of the form u(r,t) = 9 (r/t) are 
well-known to exist, as a consequence of the invariance of (|2 . 16[) under the scaling 
(r, t) 1 — !• (Ar, \t) , A e E. More precisely, from (|2.16p 8 clearly has to satisfy 

ce' (c)-(A(6(c)))' = o 

and this yields, formally, to 

u(r,t) = e(r/t) = (A)- 1 (r/t), 

as long as, (A ) _1 is well defined at least in a sufficiently small real interval. 
In the former case, one can search for particular solutions of (|2 . 1 6[) in the form 
u (r, t) = rj (r — ct) , where c € M is the speed of the travelling wave and rj has to 
be determined. Substituting the expression for u(r,t) in (|2.16p . we deduce that c 
is an eigenvalue (with 77 as eigenvector) for 

(— c + a (rj (r — ct))) rj (r — ct) = 0. 
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Under appropriate structural conditions on a (•) (see, e.g., Section 3), this eigenvalue 
problem is strictly hyperbolic (cf., e.g., [HI Chapter 11]) with speed c = c{rj) > 
0, hence a wave-like solution rj = r){r — ct) to (|2.16|) can be found. This is a 
unidirectional wave which travels into the region O. We can now map back to our 
original coordinate system to find that u (r, t) = r\ (x — ctn) is a solution to (|2.15[) . 
In plain physical terms, the mass of fluid is carried by the wave rj into an infinitesimal 
layer near the boundary T. This wave will cease to exist after some small time since 
once inside £1, the primary process is governed by nonlinear diffusion in the porous 
medium equation (|2.10[) . 

It is easy to extend our derivation to systems of the form (|l.lj) - (jl.4p and to give 
them phyiscal interpretations. For instance, these systems also occur in the pharma- 
ceutical industry by mathematical modells for the development of blood coagulation 
treatments with specific coagulation factors [H } 128 } 132] . The systems (|l.l|) - (ll.4p are 
also motivated by diffusion processes on metric graphs and ramified spaces, which 
yield interface problems for quantum graphs with coupled dynamic boundary con- 
ditions at the nodes (see, e.g., [33] and references therein). On the other hand, 
the reaction-diffusion equations (|l.ip (for A t (ui) — diUi, i € N m ) arise as models 
for the densities u i: i € N m of substances or organisms that disperse through space 
by Brownian motion, random walks, hydronamic turbulance or similar mechanisms. 
These equations are widely used as models to account for spatial effects in ecological 
enviroments [5]- For equations (ll.lj) . the Dirichlet boundary condition (|2.1ip speci- 
fies the density Uj of species at the boundary T, with an interpretation that anything 
that reaches the boundary V of fl leaves and does not return. If uq = 0, then (12. lip 
may be interpreted as if the species suffers extiction if say the patch Ti where the 
individuals live is toxic. The (homogeneous, H = 0) Neumann boundary condition 
(|2.12p says that nothing can cross the boundary of fi. Another relation in ecological 
models is the Robin boundary condition (|2.13p with j (x) = and C — Ci £ I*, 
i S N m which can be interpreted as saying that when organisms reach the boundary 
some leave it but some do not depending on the sign of Cj. Finally, the other not so 
common condition is the Wentzell-type (dynamic) boundary condition ([2. 14|) which 
states that change in the density of individuals at Ti is a function of their flux in 
the normal direction across Ti and some other function of density if no dispersive 
effects along Ti are taken into account. Following our reasonning above, this type 
of boundary condition (|1.3p can be interpreted as saying that some individuals may 
choose to live on the patch Ti but some may not and can choose to return to the 
region Q, where spatial diffusion coupled with reaction in the bulk Q is the main 
mechanism for population movement and interaction. Suppose, for instance, that 
certain critical resources for a specific population Ui, i 6 N m , are available only on 
T±. Then m must obey the rule 

(2.17) d t Ui + diVui ■ n + h(x)ui — 0, on Ti x (0, oo) , 

which says that the density Ui diffuses (in an infinitesimal layer near Ti) toward the 
patch Ti in the direction of normal flux. Again the main mechanism for this behav- 
ior here is the influence of external forces on T on a particular population Uj. Of 
course, in this context the function h (x) plays the role of a resource density function 
on Ti, and it can generally depend also on time. In fact, it is not hard to imagine a 
typical scenario where predatory individuals are preferentially concentrated around 
valued resources on Ti where the likelihood of prey is greatest. Hence, in the more 
general case of (|2.14p the state densities ui may be also allowed to carry mass on 
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Ti in contrast to the usual Robin condition for which the mass is always zero. This 
general description (|1.3|) along the patch Ti can have substantial consequences on 
the dynamics of various ecological enviroments modelled by reaction-diffusion sys- 
tems. We give a short reasonning for this behavior as follows. In the case of a 
scalar (non-degenerate) diffusion equation (m = 1, ai (•) = di), we have shown in 
[22l [23] (say, m dimension N > 3) that problem (11.7l) - (|1.8j) posseses a finite dimen- 
sional global attractor Ad yn whose dimension is essentially of different order than 
the dimension of the global attractor _4d_n_r for the same parabolic problem (|1.7I) 
with a Dirichlet /Neumann- Robin boundary condition (|2.11j) - ([2.13[) (cf., also [2"3"]V 
In particular, the correct asymptotics for the Hausdorff and, respectively, the fractal 
dimensions of Aiyn are 

|r| |ri 

(2.18) dim g Aiyn ~ C(/i,ffi) ; ,,jv-i > dim F Aa y n ~ C (fi,gi) l ^ N _ 1 , 

{vb) {vb) 

as long as v — > + . Here, C = C (/i, g\) is a positive constant that is independent of 
the size of f2, but depends only on /i and g\, and |T| denotes the natural Lebesgue 
surface measure of T C R 2 . Note that the asymptotics for the dimension of >4d-n-r is 
actually C (/i) |f2| / {ub) N ^ 2 , as v — > + (sec, e.g., |5J suggesting that the dynamics 
on .4dyn is qualitatively different than that on _4d-n-r even though both systems 
are gradient like |23) (i.e., both problems possess a global Lyapunov function). The 
asymptotic estimates in (I2.18|) are essentially determined by the instability indices 
of a properly chosen family of (hyperbolic) equilibria u* (see, e.g., [5], 22 ). One 
achieves a lower bound like (|2.18[) by computing the dimension of the unstable 
eigenspace E u of the linearization of (|1.7[) - (|1.8p around a constant equilibrium it*. 
In this case, the linearized system possesses at least n ~ C <7i) |T| / {vb) 1 ^^ 1 (as 
v — > + ) unstable solutions. This points out once again to the destabilizing nature 
of the dynamic boundary condition (|2.17[) even when the dynamics in the bulk Q is 
essentially strictly linear (see, Appendix). We emphasize that this kind of behavior 
cannot hold for the Dirichlet/Neumann-Robin boundary condition (I2.11j) - (l2.13p . 

3. Main results 
The natural phase-space for problems of the form (|l.l[) - (11.4p is 

X si, S2 := //i(ft)©L S2 (r) - {u= M : Ml eL Sl (!!), « 2 eL S! (r)}, 

si, S2 G [1, +oo] , endowed with norm 

(3.i) \\u\\ xsl . S2 = ( [ \ Ul dx) 1/31 + ( [ \u 2 ( X )r ds x x 



si 

n / \Jr 



if 8i, S2 € [1, oo), and 

\\U\\ X °° ■= max{||ui|| L oo (n) , ||w 2 ||z,-(r)} 
- hilU-(o) + ||«2||i~(r). 
We agree to denote by X s the space X s,s . Moreover, we have 

(3.2) X s = L s (O, dp), se [1, +oo] , 

where the measure d/i = dx\Q dS x \r on f2 is defined for any measurable set yl C O 

by 

(3.3) ^(A) = \Ann\+s(AnT). 
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Identifying each function 8 G C (IT) with the vector O = Q n ) , we have that C(fi) is 

a dense subspace of X s for every s G [1, oo) and a closed subspace of X°°. In general, 
any vector 9 G X s will be of the form (^) with d x G L s (fl, dx) and 6 2 G L s (T, dS) , 
and there need not be any connection between 9\ and #2- 
Next, we set 

X S1 > S2 :=TT L Sl (fi)xTT X S1 ^ 2 , forany Sl ,s 2 G [l.+oo], 
where, for any given set X, 

JJ ^X := X x . . . x X . 

\I\ times 

The norm in the space X S1,S2 , for any si, s 2 G [1, +oo) is 

(3.4) Hl^ll^,., := X) i6/m INIi'i(n) +Z! ie j m (iKIIi'ifli) +** IKHz,*2(r)) > 

where 1? = (ui, u m ), while the norm in the space A" 00 := X 00 ' 00 is naturally 
given by 

Vit\\x<=° : = m ^{ m ax|l Ui lli~(n)'? n a x ll u illL-(n}- 

If si = S2 = s, we will simply write X s instead of X Sl - S2 . Finally, without further 
abuse of notation, we will also refer to X , "r^ = (ri , ... , r m ) , as the following Banach 
space 

endowed with the natural norm in (|3.4p . 

Let us now state our main hypotheses on the source terms /i, <?i, ft-i and nonlinear 
diffusions dj, for each i G N m . 

Conditions on <Zj : The Caratheodory functions a.; (with values in R) satisfy the 
condition: 3a j > 0, Vs^ G K such that 

(3.5) a* (sj) > Q!j |si| Pl , i G N m , 
for some nonnegative pi . 

Conditions on /j, g^, ft,, : The Caratheodory functions gi, hi (with values in R) 
satisfy the conditions: 3C/ 4 , C gi , Ch i > 0, for almost all (x, t) , Vsi G R, such that 

T, iG j m 9i (x,t,s 1 ,...,s m )s i > -J2iej m C 9i N 2 ~ 
52iel m hi ( a; '*' s i' s ™) s * ^ - Sie/ m N 2 ~ c h, 

for some nonnegative Cf,C g and C/j. 

The question of global existence and the L p -L°° smoothing property for solutions 
of p.l[) ~ p.4p can be stated for the function u = (ui, ...,u m ), as follows. We say 
that the parabolic system (|l.ll) - (|1.4p satishes Property P(ri,ra), for some finite 
Ti,r2 > 1, if, for all i G N m , any of the following conditions holds: 

(i) There exists a positive function Q, independent of initial data, and a positive 
constant 77, such that 

(3.7) sup \\t(t)\\ xrur2 <Q( V ). 

t>T]>0 

(ii) There exists a positive constant C, independent of initial data, such that 

(3.8) limsup\\lt(t)\\ xri ,r 2 <C. 

t— >oo 
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(iii) If — (wio> •••) Umo) £ i G N m , there exists a positive function Q, 

independent of initial data, such that 

(3.9) sup H^WH^ <Q enroll*-)- 

t>o 

The first result concerning the L°°-estimate for solutions of (|l.lj) - (|1.4p in the 
non-degenerate case, shows that if either one of the properties for P(si,S2) above 
holds apriori for some finite Si, S2 > 1, then it also holds for si = S2 = oo. 

Theorem 3.1. Let the assumptions h3.5\) , h3. 6]) be satisfied such that p t = 0, for 
all i G N m (i.e., il.5\) holds). Suppose that the system satisfies the 

property P(l,l)-(i) (respectively, (ii) or (iii)), then it also satisfies P(oo,oo)-(i) 
(respectively, (ii) or (iii)). 

Remark 3.1. Note that we can also consider the more general case in which ai (ui) 
is replaced by (x, i, ~it) , i G N m , that is, the equations (jl.ip are strongly coupled 
in their diffusions. In this case, we must replace assumption (|3.5[) by 

(3.10) a i (x,t,~s^)>a i \'t\ Pi , V^el'", 

for almost all {x,t), and notice that all the computations performed in the proof of 
Theorem 13. II hold automatically since |^| Pl > |s;| Pl for any s = (si, ...,s m ) G R m 
(of course, in that case pi = by assumption, for all i G N m ). Theorem 13. ll can be 
also extended to systems with p-Laplacian diffusions, i.e., 



for some g i >2,i£ N m , by following, for instance, [23] ■ 

The second result is concerned with the full degenerate case (ll.ll) - p.4[) when 
the assumptions of Theorem 13. II do not hold (in particular, if it happens that pi > 
for some i G N m ). The proof is based on a truncation technique which was 
originally developed by DeGiorgi to study the regularity of solutions to elliptic 
equations, and then extensively used by many authors to study weak solutions to 
degenerate parabolic systems, subject to the usual static boundary conditions (see, 
e.g., jTO] [27]). Here, we extend DeGiorgi's method to problems of the form (|1.1[) - 
(|1.4[) . In order to avoid additional technicalities due to the different conditions that 
one can assign on the boundary T for each Ui, i G N m , we shall focus our attention 
to the case J m = N m only (i.e., we will assume that I m — 0). In this case, we 
require that the following growth assumptions hold: 



(3.11) 




\fi(x,t, Sl ,..., Sm )\<c f (y; 

\gi (x,t,si,...,s m )\ < C g (Ei £ N„ 
for some 9i,f3 i >0 and some positive constants Cf,C g . 
Theorem 3.2. Let i3.11\) hold, and assume that 3cn,o~i > such that 

(3.12) ai|^r <ai(x,t,~s>) < ct, |^| Pl , i G N m , 
for any s = (si, s m ) € M. m . Let 

(3.13) <5:= max {2,^ + l,^ + l}, 7 := max {2, ft + 1, ^ + 1* 
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Suppose that the system U. 1\) - (L~4\ ) satisfies the property P(S, r y)-(i) (respectively, 
(ii) or (Hi)), then it also satisfies P(po,oo)-(i) (respectively, (ii) or (Hi)). In par- 
ticular, for every i G N m and T, r > such that T — 2t > 0, the following estimate 
holds: 

(3.14) sup \ui(x,t)\ 

(x,t)£[T~T,T]xTi 

<Q(l + \\^\\lS([T-2t,T]xU) + ll^lliT([T-2T,T]xr)) ' 

for some positive function Q which is independent of it, time and the initial data. 
The function Q can be computed explicitly in terms of the physical paramaters of 
the problem. 

We will now show how to deduce the property P(si, S2) , for some finite s\, S2 > 1, 
for the problem (|l.ip - (|1.2p subject to a dynamic boundary condition of the form 
(|1.4|) . We shall first consider a special case. Let T c K* -1 consists of two disjoint 
open subsets Ti and r 2 , each Ti\Ti is a 5-null subset of T and T = Ti U T 2 , such 
that Ui, i G N m satisfy (|1.3j) on Ti x (0, 00) , and 

(3.15) u t = 0, on T 2 x (0, 00) , for i G N m . 

We assume that r 2 is a set of positive surface measure, and that the nonlinearities 
fi,gi satisfy the following special form of p.6[) , that is, 



(3.16) 



X/iGN„, 9i (pl~t) s li s m) Si\si\ > — X)ieN,„ I S i I — @g 1 



for some rrii > 1, and some positive constants Cf i ,C 3i and Cj,C g > 0. Note that 
(I3.6P is equivalent to (|3.16[) for m, = 2. 

The second main result gives a A^-dissipative estimate for solutions of (II. lj) . 
(|l~3l) . (O, (l3~T5l) . 

Theorem 3.3. Suppose T 2 is a sef of positive surface measure. Let the assumptions 
\3.5\) . \3.1b}) be satisfied, and let pi > for all i G N m . Then, the system 
\ 1.3]) . Ji.^[ ), \3.15\l satisfies property P(r)-(i) for = (n, ...,r m ) with ri = rrii, 
i G N TO . Moreover, if ~vfo G A 1 *, £/ien there also exists a positive function Q, 
independent of initial data and time, such that 

sup \\t (t)\\ x ^ <Q(||^o||^). 

t>Q 



Remark 3.2. Theorem 13.31 only holds if T ^ Ti, i.e., when the boundary T 2 has 
positive measure (cf. also Remark 14.21 below) . We require different arguments for 
the case when T 2 = 0. 

We shall now derive another dissipative estimate for solutions of (|l.lj) - (ll.4[) in 
«Y~^-norm which also covers Dirichlet boundary conditions (for i G /,„.) and applies 
to the case when T 2 = 0, without enforcing any further sign restrictions on all the 
PiS (compare with the assumptions in Theorem 13. 31) . Analogous to (I3.16[) , we shall 
assume that the functions fi,gi satisfy 



(3.17) 



X/ieN m /» ( X > t> s l) •••! s m, Ci) s i \ s i\ — X^eN„, ^ Si I S i ' 

SieN m 9i ( x i Si, S m ) Si \si\ > — X^gja C a - \Sj\ P — Cg, 
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for some rBj > 1, and some real constants Cf i ,C gi and Cf,C g > 0. Moreover, 
consider the (self-adjoint) eigenvalue problem for so-called Wentzell Laplacians Aw,i 
(see [22l Appendix]), as follows: 

(3.18) — a i Aip i — Cf i ip i — Aiip in fl, 

where 

2 



at (mi - 1) 

• mi + p t , 

with a boundary condition that depends on the eigenvalue Aj explicitly, 

(3.19) a i dn<Pi — Cg i <p i =A i <p i on.r 1 , 
such that 

(3.20) <fi = on T 2 . 

Here I^ is assumed to be a set of nonnegative surface measure (the case T 2 = 
may be also allowed). Our second condition on the nonlinearities fa, gi is concerned 
with some sign assumptions on Cf i and C Bi . In particular, we assume that 

(3.21) Ai := inf A x t > 0. 

Observe that, if Cj t < 0, C 9i < 0, for all i G N TO , then the first eigenvalue h.\,i 
of (|3.18p - (l3.20p is always positive, for any i G N m . Otherwise, we note that even 
when at least one of Cf i or C 9i is positive, it may still happen that (|3 . 2 1 1) holds. In 
this sense, our system (|l.l[) - (jl,4[) becomes dissipative even when at least one of the 
terms fa,gt has the wrong sign at infinity, but the bad sign is compensated by the 
other term (see, also [31] for similar assumptions). 

Theorem 3.4. Suppose that T 2 is a set of nonnegative surface measure. Let \3. 5\) . 
3.17\j , H3. 21)) hold, and assume that pi > for all i e N m . Then, the system 
1.4 ), i3.15\) satisfies property P(^r)-(i), for r — (7*1, r m ) with ri — mi, i S N m . 

On the other hand, if Pi = for some i G N m , and if IZq G X r , then this system 

satisfies P(fr)-(ii) instead. 

Remark 3.3. All the above results can be extended to models which also incor- 
porate advection effects in the domain Q and on the boundary T (i.e., the reaction 
terms fi,g% and hi may also depend on Vu^ and VrUj, respectively). Indeed, by 
making similar assumptions to (|3.6p - (13.17l) . any integral over |Vuj| may be ab- 
sorbed by diffusion in the bulk f2 with an appropiate application of suitable Young 
and Sobolev inequalities (cf. Section 4). We will return to these questions elsewhere. 



4. Proof of main results 

4.1. Proof of Theorem [37TJ In order to justify our computations, we need to con- 
struct an approximation scheme which relies on the existence of classical (smooth) 
solutions to the non-degenerate analogue of (|l.ip - (ll.4l) (if at least one Pi ^0). One 
of the advantages of this construction is that now every (possibly, very weak) so- 
lution can be approximated by regular ones and the justification of our estimates 
for such solutions is immediate. To this end, for each e > 0, let us consider the 
following non-degenerate parabolic system: 

(4.1) d t u t - div {at ( Ui ) Vu t ) + fa (x, t, l!) = 0, in ft x (0, 00) , 
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for i = l, m, where af (s$) := ai {si + e) > O (e) > 0, subject to the following set 
of boundary conditions 

(4.2) d n Ui + hi (x, t, it) = 0, on T x (0, oo) , i G I m 
and 

(4.3) Sid t Ui + a e (ui) d n (it,-) + (x, i, u ) = 0, on T x (0, oo) , i e J m . 
We equip the system (|l.ip - (|1.3l) with the initial conditions 

(4.4) «|| t =o = u m in fi > u i|t=o = u to on r > 
for uf (0) := itf e C°° (O) , i e N m , such that 

uf (0) -> u, in i s (£1) , u\ (0) |r -> w i0 in L s (r) , 

for some given s > 1. Then, the approximate problem (|4. 1[) - (|4.4[) admits a unique 
(smooth) classical solution with 

for some i* > and each e > (see [HI [TT1 [8j [31] ) . Being pedants, we cannot apply 
the main results of [16l [31] (cf. also [17]) directly to equations (|4.1|) - (|4.4|) since the 
functions af, /», and /ij are not smooth enough. Moreover, the solutions con- 
structed this way may only exists locally in time for some interval [0, t*]. However, 
by approximating the functions af, fi, gi, hi by smooth ones, say, in C°° (M,R), 
we may apply Remark 14.11 below for the solutions of the approximate equations 
(|4.ip - (|4.4p . and deduce the existence of a globally well-defined solution on [0, T], for 
any T > 0. Indeed, an apriori global bound in Holder-norm for lt e guarantees the 
global existence of classical solutions (see, e.g., [31]). Nevertheless, even when these 
bounds are not available apriori, we may still choose to work with locally-defined 
in [0,t»] , smooth solutions ~it e , that are globally defined on R + in the lower-order 
L p -norms; this turns out to be sufficient for our purposes. Indeed, if the solution it e 
is globally-defined on [0, T] in A^-norm, then it will also be global in A' 00 -norm by 
(iii). As we shall see in this section, a global bound in A^-norm can be established 
for the same assumptions p.5[) - (l3.17[) on the nonlinearities. 

We begin with the proof of Theorem 13.11 by following similar arguments to 
[TT1 H3] for the system (jl.ip with static boundary conditions. From now on, c 
will denote a positive constant that is independent of t, e, n, it and initial data, 
which only depends on the other structural parameters of the problem. Such a 
constant may vary even from line to line. Moreover, we shall denote by Q T (n) a 
monotone nondecreasing function in n of order r, for some nonnegative constant r, 
independent of n. More precisely, Q T (n) ~ cn T as n — > +oo. We begin by showing 
that the ^"-norm of it = lt e satisfies a local recursive relation which can be used 
to perform an iterative argument. We divide the proof of Theorem 13.11 into several 
steps. 
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Step 1 (The basic energy estimate in X n+1 ). We multiply (|4.1I) by \ui\ n ^ Ui, 
n > 1, and integrate over f2, for each i G N m . We obtain 

(4.5) | Il«ill2&(n) + (/< (*. f > ^ ) . kl"" 1 ^) i2(n) 

+ n J al (ui) \Vui\ 2 \ui\ n 1 dx 
n 

Similarly, we multiply (|4.2j) and (|4.3|) by |itj| n 1 itj and integrate each relation over 
T. We have 

(4.6) (n+~l) ~L ll^ll^+Kr) + y a i i u i) d ^ u i \ui\ n ^ ndS 

\ / x, (r) 

= 0, 

for each i G J m , and 
(4.7) 

9 n Mi l^iT" u «d5 + (a- /ij (a;,i, it) , |«i| n- Uij ^ = 0, i G J m . 
\ / i 2 (r) 

In the case when (|1.2[) is replaced by a Dirichlet boundary condition for Uj, i G Z m , 
equation (|4.7p still holds since /ij = in that case. 

Let us first observe that, in light of assumption (|3.5[) . we have a\ (sj) > cti, 
Vs; G R, e > (recall that y>j = 0), which immediately implies 

(4.8) f af («i) |Vw,| 2 luil"" 1 cfe > a, / |V^| 2 \u l \ n ~ 1 dx, i G N m . 

Moreover, on account of of the assumptions (|3.6[) for /j , and a basic application 
of Holder and Young inequalities, wc deduce 

and 

(4.10) Xl ie , 7m 0>*>^) > l Ui l"' lMi ) L 2 (r) - ~ c 5I ieJm 5i H^ll^+Mr) _c - 

In order to estimate the source terms involving hi on the boundary (|4.7p . we 
need the following lemma which allows us to control surface integrals in terms of 
volume integrals (see, e.g., [23 J. 

Lemma 4.1. Let n > 1, p > 0, s > — 1. T/ien /or every e > 0, £/iere ZioWs 

C 



M-«<.(. +B ) J [ i |v^|^-* + 7 (. + ,)(Mja w) + 1 ) 

/or some positive constant C = C (p, s) independent of u, e and n. 

Applying Lemma |4. II to u — Ui, with s = 1 and p = 0, i 6 J m , we have 

(4.11) / |«i| 1+T *d5< e(n + l) / |V%| 2 |u|" _1 & 

Jr Jn 

+ j(n + l)(|K||S+i 1(n) + l). 
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Thus, in light of the third assumption of (13.61) . using (|3.5p . we find 

> -E ieJro \ a i c hi Il«ill2"+i(r) IK-II2»-i(d) 

> -c^ ieJ Q * Il«»ll2»+Hr) ~ c 
which can be bounded below by 

(4.13) ~ C£ J2 £I a d n + l ) J |Vw*| 2 lu^dx 

-iE <e/m ( n + 1 )(ii«*n^(n)+ 1 )' 

exploiting the estimate (|4.11[) . Choose now e = e„ > in (|4.13p such that 

(4.14) eri = m ax-f^-, Vn>l, 

«eN m 2 (n + 1) 

and note that e„ < c, uniformly as n — > oo. Summing the equations (|4.5[) (respec- 
tively, (|4.6I) and (|4.7|l ) over the sets i G N m (respectively, i G J m and i G I m ), 
then adding the relations that we obtain, on account of (|4.9 |) - (|4.10|) . (|4.12| - (|4.13p 
we deduce 



(4.15) ^\&+cn(n+l)^ i&m J n \V 
Qa(n)(l|t?ll*tix 



Ui\ \uj\ n 1 etc 



< g 2 (n) ^ii^h^ + 1 

for any n > 1. Here, the function Q2 (n) ~ as n — > 00. 
Step 2 (The local relation). Set n k = 2 k — 1, k > 0, and define 

(4-16) y*(t):=||^WKlx, 
for all k > 0. Then, using the basic identity for u = Ui, 

1 



(4.17) / |Vu| 2 M"" 1 dx = 





n + l 


(n + l) X 


V|u| 2 





from (|4.15j) it holds 

(4.18) j t y k (t) + £. gNm 7nfc ^ I v H^fd* < Q 2 (n k ) (y k (t) + 1) , 

for all k > 0, where < 7 < 7„ fc ~ c, as rift — > 00. Let be two positive 
constants such that t — fi/n k > 0. Their precise values will be chosen later. We 
claim that 

(4.19) y k (t) < M k (t, ft) := co (ji) (n k ) a ( sup y k -i («) + l) 9fc , Vfc > 1, 

s>t—fi/nic 

where Co, <r are positive constants independent of k, and k > 1 is a bounded 
sequence for all fc. The constant Co (/i) is bounded if /1 is bounded away from zero. 

We will now prove (|4.19p when 2 < N. The case N < 2 requires only minor 
adjustments. We will follow an argument similar to the proof of [HI Theorem 2.3] 
(cf. also [23]). For each k > 0, we define 

_ JV(n fc + l)-(JV-2)(l + n fc ) _ _ 

n ' k - N(n k + l)-(N-2)(l + n k . 1 )' Sl ' k - ^ 
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We aim to estimate the term on the right-hand side of (|4. 1 5[) in terms of the 
^ 1+n *- 1 -norm of u. First, Holder and Sobolev inequalities (with the equivalent 
norm of Sobolev spaces in W 12 (O) C i Pa (fl), p s = 2N/ (N - 2)) yield 



(4.20) 



\uA 1+nk dx< 



N - 2 



< C 



,M 



|u,l 1+nfc - 1 d a ; 



r it k 



\uA 1+nk - x dx 



\u,\ 1+nk dx 



with s.^fc := Si_kN/ (N — 2) e (0, 1). Applying Young's inequality on the right-hand 

side of (j4.20j) . we get 

(4.21) 



■it; 



1+rifc 



dx < 



dx + Q T1 (n k ) / \u t \ 



il+n fc _i 



dx 



for some positive constant t± independent of rik, and where 

zi,fc := r it k/ (1 - Si,fc) > 1 

is bounded for all k. Note that we can choose to be some fixed positive number 
since Q T2 also depends on 7 ~ c. 

To treat the boundary terms on the right-hand side of (|4.15|) . we define for k > 0, 

(N - 1) (n fc + 1) - (JV - 2) (1 + n k ) 



Ui.k 



-, %i,k := 1 - Vi.k- 



(N - 1) (n* + 1) - (iV - 2) (1 + n k -i) 

On account of Holder and Sobolev inequalities (e.g., W 1 - 2 (Q) C L qs (T) , q s 
2 (N - 1) / (N - 2)), we obtain 



(4.22) 



|«<| 1+n * dS<c 



< c 



(jV-l)(i fc + l) 



Vliii 



lUil 1 ^" 1 dS 



dS 



dx 



Vi.k 



\u i \ 1+n "- 1 dS 



\uA 1+nk dx 



with Wi t k ■— %i,k (N — 1) / (N — 2). Since x^k E (0, 1), we can apply Young's in- 
equality on the right-hand side of (|4.22[) . use the estimate for the L 1+Tlk (fi)-norm 
of m from (I4.20| in order to deduce the following estimate: 

li.k 



(4.23) 



,1+n, 



dS<^ 

4 



V|u ? ; 



dx+Q T2 (rife) / \ui\ 



il+n fc _i 



dx 



for some positive constant T2 depending on t\, but which is independent of n k , and 
where 

H,k : = T\ = — V - 1 

is bounded for all k > 0. Inserting estimates (|4.2ip - (|4.23[) on the right-hand side of 
(|4.18l) . we obtain the following inequality: 

2 



(4.24) d t y k (t) + 



ieK, 



V|u, 



dx<c(n k ) ai + !)"*, 
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where c, o\ are positive constants independent of fc, and 

9k ■= maxfmaxfz^} ,max{Z i]fe }) > 1 

i i 

is a bounded sequence for all k. 

We are now ready to prove (|4.19[) using (14.241) . To this end, let C, (s) be a positive 
function C : R + -> [0,1] such that £ (s) = for s 6 [0,t-/V«fe]> C (s) = 1 if 
s e [t, +oo) and \d£/ds\ < n k /^i, if s G (i — M/ n fci We define Z fc (s) = £ (s) 34 (s) 
and notice that 

Zk ( s ) < -^y k ( s ) + c s — 34 s . 

as /i as 

Combining this estimate with (|4.18|) . (|4.2ip . (|4.23l) and noticing that Z k < 34 > we 
deduce the following estimate for Z k : 

(4.25) —Z k (a) + C (/x) (s) < M fc (i, /*) , for all sg[t- /i/n fc , +co) , 
as 

for some positive constant C independent of k. Integrating (|4.25|) with respect to 
s from t — n/nk to t and taking into account the fact that Z k (t — fi/n k ) = 0, we 
obtain that 

y k {t) = z k (t) < M k (t, fi) (i - e - c ") , 

which proves the claim (|4.19j) . 

Step 3 (The iterative argument). Let now t > r > be given; define fi = 2(r — r), 
to = t and ife = t^-i — l^/n k , k > 1. Using ()4.19|) . we have 

(4.26) sup 34 (t) < c {n k ) a (sup 34-1 (s) + 1)<\ fc > 1. 

t>tfc-l s>t fc 

Here Co = Co (/i) depends only on /i. Now let us define 

(4.27) C:= sup (y (s) + l)= sup (||^ (s)||^ + 1) . 

s>ti— r s>ti— r 

Thus, we can iterate in (|4.26l) with respect to k > 1 and obtain that 
(4.28) 

sup y k (t) < (co<) (conU) 8 * (conZ^Y* 8 *- 1 ■ ... • (cdO**'"- 1 -* (C)^ 
t>t)=-i 

<<£w B >> (c) ik , 

where £ fc := 6k6k-i—6o, and 

(4.29) A fc := 1 + fc + fc fc _i + ... + 9 k d k -x...6 , 

(4.30) B k :=k + 9 k (k - 1) + fe fc _! (fc - 2) + ... + Mfc-i-Ao- 
We can easily show that 

x 1 x ' 

(4.31) A k < (ci + n fe ) V ■ and B k < (c 2 + n k ) V 

^— ' Ci +71- ^-^ 



i=i 



i " C2 + n J 



for some positive constants ci,C2 independent of fc, /i. Therefore, since 
(4.32) su V y k {t)< sup 34(*) ^c^" 5 * (Cf k 

t>t *>**-! 
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and the series in (|4.31[) are convergent, we can take the 1 + n^-root on both sides 
of (|4T32]) and let k ->• +00. We deduce 

sup \\lt (t)\\ XOB < lim sup (y k (t)) 1 /(i+»») 5 

t>t =r' ^+°°f>*o 

which, on account of (|4.32[) . yields 

(4.33) sup \\!t(t)\\ X - <C( M )(C) 1/C3 , 

t>t =r' 

for some positive constant C3 independent of t, k, it , e, /i, and initial data. Note 
that C depends on r (see Q4.27|l ). 

Step 4 (The final argument). Let us first assume that Property P(l)-(i) holds. 
Then we already know that the A^-norm of it (t) is bounded independently of the 
initial data, for each t > r. Therefore, from (|4.33j) we also obtain the claim for 
the A'^-norm of it (t), i.e., Property P(oo)-(i) holds as well. If, on the other 
hand, Property P(l)-(ii) holds, we can choose r = r + 2/z with fj, = 1 so that C 
and C (/Lt) are bounded uniformly with respect to initial data as r — > 00. Hence, 
Property P(oo)-(ii) is also satisfied by letting r — > 00 in (|4.33j) . In order to show 
the final property (iii), taking advantage of the fact that the initial data ito S X 00 , 
it suffices to note that in place of the inequality (I4.19[) . we may use instead the 
inequality 

y k (t) <Q(\\lt \\ xoa , sup M k {t,n)), 

which is an immediate consequence of (|4.24j) . Arguing analogously as in |31[ Lemma 
5.5.3], we obtain the claim. The proof of Theorem 14.331 is now complete. 



Remark 4.1. It was proven in [3TJ Section 5] that maximal L p -regularity for (|4.1jl - 
(|4.4|) can be used to reduce the question of global existence of the solutions Tt e in a 
space of maximal regularity, to the boundedness of lt e in a Holder norm C 0,13 (f2) , 
/? > 0. It should be possible to prove, under the natural assumptions of Theorem 
13.11 that every classical solution of problem (|4.1I) - (I4.4[) is globally Holder continuous 
on Q. Establishing global Holder continuity for solutions to systems with dynamic 
boundary conditions requires a more detailed analysis, involving careful local esti- 
mates of the solution near the boundary. Of course, as in the case of Dirichlet/Robin 
boundary conditions for (|4. 1 [) (see, e.g., [12]), these Holder bounds should apriori 
depend on the L°°-norm of the solution. Thus, our analysis constitutes only the 
first step in proving boundedness in Holder norm C '^ (f2) . This question remains 
open for now. 

4.2. Proof of Theorem 13.21 We shall divide the proof into several steps. As 
in Section 4.1, we can justify our computations by exploiting the approximation 
scheme (|4.1|) - (|4.4|) . As before, c will denote a positive constant that is independent 
of t, e, n, it and initial data, which only depends on the other structural parameters 
of the problem. Such a constant may vary even from line to line. Without loss of 
generality, we may assume that 5^ = 1, for all i € J m = N m . 

Let T,t and L be positive numbers such that T — 2t > and L > 1. We set 
to = T — 2t and define the sequences 

t n = t n -i + k n =L(2--), for all n > 1. 
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Consider the (smooth) cut-off functions r\ n G C 1 (R, [0, 1]) with the property that 

= i 0, t <tn _l. 

Next, denote Q« := / n x O, where /„ := [i n _i,T], and the sets 

:= {fa, *) S I n x : (x, t) > k n ) , A r l n := {(x, t) e I n x T : u, (x, t) > k n } . 
Let Ai tn = Af n U Af n , and note that 

Ai : n — \\X, t) € Q n '. Ui (x, t) > k n } . 

Finally, we denote by \A^ n \ the (N + 1-dimensional) Lebesgue measure of the set 
Af n , and by \A^ n , the (TV-dimensional) Lebesgue measure of Af n , respectively. 
We note that, according to (|3.2j) - (|3.3|) . we have 



|Q n | = |/ n | M (n) = |/„|(|fi| + |r|), 

and we can do so similarly for the set A i n . 

Step 1. (The energy inequality). We define the truncated functions 

u i>n (x,t) := max{ui (x,t) - k n ,0} = (u, - k n ) + . 

We begin by multiplying equation (|4.1[) by Ui tn r] n (t) and integrating the result- 
ing identity over /„ x Cl. Then, we multiply (|4.3[) by Uj )ra ?7„ (t) and integrate over 
I n x F. Adding as usual (cf, e.g., Section 4.1), then exploiting the growth assump- 
tions (|3.1ip on fi and Qi, the fact that rj n (t) < 2"/t, we obtain after standard 
transformations 
(4.34) 

2 *\ J_ i /„,2 L i // „ . f„. .\ IV7 t„,. „ \ rtl 2 



max / u i n (x, t)dx + / u ln (x, t) dS + a% (ui) |V {u % , n rj n ) {x, t)\ dxdt 

\Jn ' Jr ) JJi n xn 

' ' i I (-^r u ln (x, t) r} n + ^2 jeN \uj (x, t)f u hn (x, t) i] 2 n (t) + Ui, n (x, t) r£j dxdt 

(on \ 
— u ln ( x > *) Vn + X^- eN \ u i ( x ^)f^,n (x,t) i}\ + Ui, n (x,t)rfcj dSdt, 

where we have set 9 := maXj e N m Qi and /3 :— maxj e N m f3 i . Now, we wish to estimate 
the terms on the right-hand side of ()4.34j) . To this end, set 

A n := (J k >iA% n , jF n := U fe >i^ n 

and note that on A^\Af n and A„\Af n , respectively, we have m (x,t) < k n < 2L 
and Ui (x,t), r < k n < 2L, respectively. Therefore, 

(4.35) / / dxdt < cL e+1 [ [ (1) dxdt < cL e+1 V \A?\ , 

and, analogously, for the trace of it; we have 
(4.36) 

/ / \u,f +1 dSdt < cL f3+1 1 1 (1) dSdt < cL p+1 J2 jeN K«| , 
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Since k n ~ L, as n — > oo, it is easy to see that the following inequalities hold: 
f \A% n \ < ck« +1 \A%\ < cjf A? \u,\ a+1 dxdt, 

\ L@ +1 \Al n \ < ckt 1 < c %" N ^ dSdL 

From these estimates, we thus find that 

(4.37) JJ^ \uj\ 6 Ui , n dxdt < JJ (\u/ +1 + dxdt 

<cV [[ \uA e+1 dxdt 

and, similarly, 

(4.38) jj \ Uj f Ui, n dSdt < fj (\ujf +1 + dxdt 

Hence, using the above inequalities (|4.37l) - (|4.38l) on the right-hand side of (|4.34[) . 
and summing the resulting relation over i £ N m , we deduce 

(4.39) 

m E f* (E, eNm J u l™ *) dx + E ieNm J T dS 



lei 



2 n c /* f 2 n c f f 



where <5 and 7 are defined as in f|3 . 13|) . Here we have also used assumption (|3.5p . 
Step 2. (Additional estimates). From the definition of k n , we see that 1— fc n /fe n +i > 
2 -(n+2) j w hich yields 

(4.40) / / H 5 cfadi < 2^+^ f f \ Ul f (l - -^-) dxdt 

< 2 n5 c JJ ^ ( Ul -k n ) S + dxdt. 
Moreover, the same argument gives 

(4.41) / / Unp dSdt < 2" 7 c // (m - k n )\ dSdt. 

On the other hand, since on Af n+1 UAf n+1 , we have (ui — k n ) + > k n+ i — k n , there 
holds 

(4.42) JJ^ ( Ui - k n ) S + dxdt > {k n+l - k n ) S \A% n+1 \ > \A? >n+1 \, 
JJ at ~ MI d Sdt > (k n+1 - k n )' ( \Al n+1 \ > \Al n+1 \ . 
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Because of these two inequalities (|4.42l) . for any positive number A such that, if 
A < S and A < 7, on account of Holder's inequality, it also holds 

(4.43) 



A n \ J J 



(m - fc„+i)+ dxdt < ( // {ui-k n+l f + dxdt\ \Af n+l \ X 



L s-x 

and 
(4.44) 



C 2"(«5-A) 

< " r r„ A / / (ui - k n ), dxdt, 



A/7 



(Ui- kn+x^dSdt < J ^ (ui - k n+1 )'ldxdtj l^n+ij 1 * h 

< L , ( _ x J J Ar - k n y + dSdt. 

Next, we will collect some useful inequalities which follow from the following 

well-known embeddings: H 1 (Q) C p s := 2N/ (N ~ 2) and H 1 (ft) C U' [T) , 

<? s := 2(N- l)/(iV-2). We give the argument for N > 2, the case AT < 2 

can be treated analogously. Suppressing the dependance on the subscript n for 

the moment, from Holder's inequality and these Sobolev embeddings, for every 
v Mi e w i,2 x we haye 

JV— 2 _2_ 

(4.45) ^ v s dxdt< J^J v mPs dx^J " (J v 2 dx^ dt 

< (JJ^ (\Vv Mi f + \v\ Mi ^dxdt \ x (max J v 2 dx 
where, for each Mj > 0, we have set 

4 

s = 2Mi H . 

1 N 

Similarly, for each M, > and Z = 2M l + 2/ (A 7 " — 1) , we have 

iV-2 1 

(4.46) ^ v s dS < J^J v^'ds)^ v 2 ds) N ~ 1 dt 

- UL ( |WV/f + H ") dxdt ) x fe x I/ ds )^ ■ 

Exploiting (|4.45p - (|4.46j) with M % = pi/2 + 1,1 = I n+1 , v = - fc„+i) + , we get 

(4.47) ff ( Ui -k n+1 ) s + dxdt 



< 



(_//" q () Ui \ Pi l^^^+il 2 + u t!n+i) dxdtj x ^max J u 2 in+1 dx 
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and 
(4.48) 



(ui - k n+ i)\ dSdt 

-{^jjj xn (\ u i\ Pi \^ u i,n+i\ 2 + "i,n+i) dxdtj x ^max J u 2 in+1 dx^j 
Finally, from (|4.39[) we see that estimates (|4.4Q|) - (|4.41[) yield the following inequality 
(4.49) max (]T ^ jf <„ (x, t) dx + ]T ^ jf uf,„ (*, t) ds' 

+ „, // a 2 : |ui| Pi |V(Mi !rl ?7„)| 2 da;dt 



< / / (ui — fc n -i)i dxdt 



A 1 



Step 3. (The iterative argument). We continue our main argument by first recalling 
the following result (see, e.g., [TU1 Lemma 4.1]). 

Lemma 4.2. Let {y n } be a sequence of positive numbers such that it satisfies 

(4.50) y n+1 < cb n y l n +K , 

for some constants C,b, k > 0. If 3^o < C^ 1, lK b~ X ' Ik , then y n — » as n -A oo. 
Define 

X,n := j^-j- («i - &n)+ dxdt + J J (ui - k n )1 dSdt^j , 

where we recall that Q n — I n x f2 and 

|Q n | = \I„ x fi| + |7 n x r| . 

Set y n — Y^ILi y%.n- The goal now is to show that the sequence {34} satisfies a 
recursive relation of the form (|4.50[) . First, using the definition of y n , we can rewrite 
(|4.49p as the following inequality: 

( 4 - 51 ) (E ieNm j n U ln+1 (X, t) dx + j v <„+! {X, t) dS^j 

+ // a^u.f 1 \V (u hn+1 n n+1 )\ 2 dxdt 

< c2 («+i)(p+i) r -i|Q ii+1 |3; n; 



M, 



for p := max (7, 8) > 1. Secondly, applying (|4.43[) to the bulk integral over u, in+1 
(where M{ :— Pi/2 + 1), which occurs in the integrals in (|4.47l) - (|4.48[) . and then 
using (|4.51l) to estimate the second terms in those products, we obtain 

(4.52) // ( Ui - k n+1 ) s , dxdt 



<c[ — + TS _ M _ ) \Q n+1 \ y n ( — — \Q n +i\y n 
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and 
(4.53) 



(ui - fc„+i), dSdt 



< c 



T 



\Qn+l I D^ra 



£7 -Mi 

Holder's inequality applied to 3^,n+i yields 
1 



c2 n P 



IQn+lIX 



(4.54) 



|Qn+l| 
1 



(ttj - k n +i), dSdt 









1 A n 1 




1 i,n+l | 








4 r 1 

■"■i,n+l \ 



1-6/s 



i- 7 /; 



|Qn+l| 

Inserting the estimates for |Af„ +1 | and |^4f„ + i|, respectively, from (|4.42[) . on the 
right-hand side of (|4.54[) . we deduce 



yi,n + i<c\Q n+1 \-yl +N ° ( — + £A -_ Mi 



(5/s 



c2 n ^ 



c|Q„ H 



2 1 
77 



f -D ! U, (JV - 1)i 1 

r L 7 ~ Mi 



7 /i 



c2 n P\ («-!)' / 2"T 



l-7/i 



Henceforth, by setting 

« := k(5,7) = 



f 25. 7 ^ ;f y > i 

if y„ < i, 



max 

mm (l^' (iV-l)i 

the above inequality yields the recursive relation 

y n +i < Cb n y l + K , with K > o, 

where C ~ L~ a depends on r -1 and b ~ 2^, for some positive constants cr,C 
depending on <5, 7, s, Z. Therefore, if we choose L > 1 sufficiently large so there 
holds 

< c^V^-i/* 2 < L°'/ k 6- 1 / k2 ) 
then by Lemma |4~21 it follows that y n — > as n — > oo. This implies that 

sup (a;, t) < lim fc„ < 2L. 

(x,t)e[T-r,T]xn ,wo ° 

In order to estimate Ui (x, t) from below it suffices to apply the result just ob- 
tained to the functions Ui (x,t) — —Ui (x,t), which satisfies a system of the same 
type as for u, (x, t) , with nonlinearities Hi (x, t, ui) = — a, (x, t, —ui) , fi (x, t, ui) — 
—fi(x,t,—ui) and gi(x,t,Ui) = —gi(x,t,—ui), respectively. These functions are 
subject to the same conditions (|3.12p . (|3 . 1 1 [) . This yields the desired estimate 
(|3.14[) . Finally, we may conclude that if T is sufficiently large, we can take r = 1 in 
(|3. 14|) , which also immediately gives the first conclusion in the theorem. The proof 
is finished. 
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4.3. Proof of Theorems 13.31 13.41 In order to justify our computations for prob- 
lem ()1.3|) . (|1.4|) . (|3.15[) . it is not clear how to use the scheme introduced at 
the beginning of the section due to the nature of the boundary domain (indeed, 
Ti n r 2 ^ in that case, so we cannot exploit maximal regularity theory to con- 
struct smooth solutions unless T 2 = 0). However, the proof can be based on the 
application of a Galerkin approximation scheme which is not standard due to the 
nature of the boundary condition (|1.3[) . We refer the reader to, e.g., [7J, for further 
details, where a system of reaction-diffusion equations for the phase-field equations 
with dynamic boundary conditions were considered (cf. also [2 6) in a degenerate 
case). 

We begin the proof of Theorem 13.31 with a result for the eigenvalue problem 
for so-called Wentzell Laplacian (see jUJ Appendix]). More precisely, let us 
consider the equation 

(4.55) - Aip = Atp in Q, 

with a boundary condition that depends on the eigenvalue A explicitly, 

(4.56) d n <p = A<p<mT 1 , 
such that 

(4.57) d n <p + ip = on Y 2 . 

(recall that r 2 is assumed to be a set of positive measure and that (|4.57j) holds 
where a Dirichlet boundary condition for u = ut is satisfied on r 2 x K+). Such 
a function <p will be called an eigenfunction associated with A and the set of all 
eigenvalues A of (|4.55[) - (I4.57|) will be denoted by Aj, j G N. Let (p x £ C 2 (H) and 
Ai, denote the principal eigenfunction and eigenvalue of (|4.55p - (|4.57j) . respectively. 
We have the following. 

Proposition 4.1. For the spectral problem |^.5<5| )- |^.57| ) ; Ai > is simple and 
(p x > in Q. 

Proof. Using the standard characterization for the eigenvalues Aj of Aw (see, e.g., 
|22j). we obtain that the following min-max principle holds: 

(4.58) Aj — min max R w Up, <p) , j G N, 

dim Yj —j 

where the Rayleigh quotient Rw, for the (boundary perturbed) Wentzell operators 
Aw, is given by 



(4.59) Rwfa<p) ■= "' T "' Z - ; 2 r ' L2(r2) , O^yGff 1 ^). 

Exploiting a well-known Friedrichs-Poincare's inequality, we have Rw (f, <f) > 
°w ll'/'llx 2 ' fo r some positive constant cw, which implies that Aj > 0, for any j G N. 
By the maximum principle, ip 1 is positive in fl since T 2 has positive surface mea- 
sure. The fact that Ai is simple, follows again from the maximum principle (see, 
e.g., S). " □ 

We are now ready to give the proof of Theorem 13.31 Without loss of generality, 
we can take Si = 1, for all i G J m . We multiply (jl.ll) by lu^™* -1 sgn (uj) ipi, and 
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integrate over fi, for each i G N m . We obtain 

(4.60) — "TT / \ui\ m ' Vidx + (fi (x,t, it) , l-Uil™ 1 " 1 sgn (ui) <p x \ 

mi at Jq \ I l 2 (Q) 

div (ai (ui) Vuj) \ui\ m ' 1 sgn (ui) ip-^dx 

o 

= 0. 

Similarly, we multiply (|1.3|) by |iti| m * 1 sgn (ui) (fi and integrate the relation over 
T (recall that (f37T5]) holds over T2). We have 



1 d 



(4.61) ~~k I \ u i\" 1 ' fidS + I a,i (m) dnUt \ui\ m ' 1 sgn(u l )^ 1 dS 
nii at i /p i Jp 

+ Igi (x,t, it) , l^l" 1 ' -1 sgn (ui) ipA 
\ I £ 2 (ri) 

= 0, 

for each i G N m . Consider the following real- valued function 

(4.62) = ^-\u t (x,t)r. 

Integrating by parts in (|4.60p . then using (14.61[) . on account of the following com- 
putation 

div (cii (ui) Vuj) |ui| m * 1 sgn (ui) f-^dx 

-(mi -I) I ai (ui) IVw^ 2 \ui\ m *~~ 2 ip-^dx 
Jn 

ai (ui) \ui\ mz 1 sgn (ui) Viti ■ Vip 1 dx 

+ J a i (ui)d n Ui\ui\ mi ~ 1 sgn(ut) Lp-^dS, 
we deduce the following inequality 

(4.63) d t _£ (it (t))tp 1 dfi + ^2- (nii - 1) / a,; (ui) \Vm\ 2 \ui\ mi 2 kp x dx 



+ ^ / a-* («;) M™ 1 1 s#n ( u %) Viti • Vip-^dx 
< c £ (it (t)) (fidfj, + c. 



Here we have employed (|3.16|) to estimate the terms involving fi,gi in (j4. 60|> - f|4.61[) . 
Let us now estimate the third integral term on the left-hand side of (14. 63[) . Exploit- 
ing the assumption Q3.5P on a;, i G N m , we deduce for u = Ui, p = Pi, M = m; 
that 

(4.64) / a(u)|u| p 1 sgn (u) Vu • Vip x dx > / |u| p+M 1 sgn (u) Vu ■ Vip^dx 
Jn Jn 

= / Va (u) ■ V<p x dx, 
Jn 



where we have set 



a(s)\s\ p - 1 ds> c\u\ M+p . 
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Integrating by parts in (|4.64|) once more and noting that a (0) = 0, we obtain 
(4.65) / Va (u) ■ Wcpidx = / a (u) dnip^S — / a(w)A(p 1 da; 



a (u) Aip^x 

a(u)Aiip 1 dS + I a{u) ki^p^x 

ri 



>Ai / \u\ M+p ip^, 
Jn 

since (Ai, ip^ satisfies (|4.55|) - (|4.57[) . Inserting the above estimates in (I4.63p , we get 
the following inequality 

(4.66) d t [e (^(t))^^ + AiV. N [ \ Ui r +Pi cp^fi 

Jn ^ !6R '" m 



< c l_£ (~u (t))tp 1 d f i + c. 
Jn 

Since all p^s are positive, we can absorb the integral term on the right-hand side 
of (|4.66p . using the Young's inequality as follows: 

y \u i \ mt < e y \u % r +p -+c e , 

for a sufficiently small e G (0,Ai) and some positive constant C e , independent of 
Ui,t. Moreover, setting v — minj e isj m (rrii/pi) + 1 > 1, we immediatelly have from 
the above inequality, that 

(4.67) / (£(lt(t))f Vl d^<cy [ Imr+^^dfx + c, 
Jn ^eN m 

for some positive constant c independent of u, t and initial data. Using (|4.67[) . we 
see that (|4.66p yields the following estimate 

(4.68) d t [_£ {it (t)) tfidp + c [_(£(lt [t)))" tprfu < c, 

Jn Jn 

by an appropriate choice of e < Ai/2. By normalizing the eigenf unction ip 1 in (|4.68l) 
such that [Ij^ifo d/i) = •"■> 011 accoun t of Jensen's inequality, it follows that 

(7 erf®) ^dfi) < he(it(t))) v tp 1 d l i, 

\Jn J Jn 

which gives the following estimate: 

(4.69) d t Y(t) + Cl (Y(t)) u <c 2 , 
for some positive constants ci,C2, where we have set 

Y(t) := / £{lt{t)) Vl dii. 



We can now use the Gronwall's inequality (see, e.g., [37l Chapter III, Lemma 5.1]), 
applied to (|4.69[) to deduce that 

(4.70) Y (t)<(^\ +(d(i/-l)i)-^, Vi>0, 
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which yields the desired claim. The proof of the theorem is complete. 

Remark 4.2. In the case when Tto e X 1 *, F(0) = lim i _>. + Y (t) is finite, 
similar argument to |371 Chapter III, Lemma 5.1] gives 

Y (t) < max (y (0) , f—\ X , V* > 0. 



Thus, the second assertion in Theorem 13.31 also follows. We also note that the above 
argument relies entirely on the fact that the boundary T2 has positive measure and 
this gives Ai > 0. The proof fails to work if for instance, F = Fj (i.e., when T2 = 0). 
We shall require different arguments for this case (see below). Finally, if at least 
one pi = 0, for some i G N m , the above argument can still be used to derive the 
following bound 

sup (t)||^ <Q(||l*oLr* e ct ). 
t>o 

Indeed, this follows from a standard application of Gronwall's inequality to (|4.66l) . 



We now continue with the proof of Theorem 13.41 As in the proof of Theorem 
3.3[ we multiply (jl.ip by \ui\ m ' 1 sgn (ui) , and integrate over fl, for each i £ N m . 
Then, we multiply both equations (|1.3[) and (jl.2p by \ui\ m ' sgn (u^ and integrate 
the relations that we obtain over T. Analogously to (|4.60l) - (|4.61l) and arguing in a 
standard way as in (|4.7[) , we deduce the following identity: 



Ui\ 1 dx + / \ui\ 1 dS 
n Jr 



(4-71) #E - 



+ E l6Nm {9i{x,t,f),\ui\ mi l sgn{ui) /L ^ 

= - ^ (rrii - 1) / a, (a;,t, "zf) |Vuj| 2 |ztj| ms_2 dx. 

By assumption (|3.5[) . we can estimate the term on the right-hand side of (|4.71[) as 
follows: 



(4.72) f aiit^^Ui^Ui^^dx^CH [ |Vu 2 | 2 \u t 
Jn Jn 



rrii-\-pi—2 



a, 



mi + pi 



dx 

2 







)7„ 


V|«| 


2 









To estimate the nonlinear terms on the left-hand side of (|4.71j) . we may exploit 
(|3~T71) . On account of (|4~72l . we have 

(4.73) |E i ^^f/w m< ^+/j^r^ 



dt ^— 'ieN m to 



a,- 



f 2 v k 1) / 


m-i+Pi 






2 


ym.+pij J n 







2 

rZ.r 



- E ieNm / o w m<+K dx + c * / r w mi+Pi dS 



< c. 
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By assumption (|3.21[) . it follows that, for all i £ N m , it holds 
(4-74) at ||V^||^ 2(0) - C fi IlipJI^fn) - C g% I^J^r) 

> A M (hiWl^n) + ll<ftll^(r)) , 

for all ip i € H 1 (Q). Thus, by choosing ip i = \u i \ lja%+v% ^ 2 in (|4.74|) . and recalling 
(|4.62p . from (|4.73l) . we obtain the following inequality: 

(4.75) d t j_£ {it (*)) + Ai J2 i&m (J N™ I+Pl dx + C m J \ui\ mi+Vi dS^j 

< c. 

Let us assume that pi > 0, for all i € N m . Arguing now as in the proof of Theorem 
13.31 (see (|4.66[) - (|4.69[1 ). it is not hard to see that we arrive at the following inequality 

(4.76) d t J_ £(lt(t))d f j, + c(J_ £(lt(t))d f A < c, 

where v = mini g N m (m-i/pi) + 1 > 1. Thus, we can use the Gronwall's inequality as 
before (see (|4.70p ) to derive the estimate 

(4.77) J_£ {it (t)) dfi < c (l + t~-^^ , Vt>0. 

Hence, the first claim of the theorem follows from (|4.77j) . On the other hand, if at 
least one pi = for some i G N m , we obtain the following analogue of ()4.76|) : 

(4.78) d t [_£ {it (t)) d f i + c [_£(Tt (t)) dfi < c, 

Jn Jn 

which yields the second claim of the theorem once more on account of Gronwall's 
inequality. In particular, there exists a positive function Q, independent of initial 
data and time, such that 

(4.79) sup ||t? (t)\\ x + < Q (ll^oll^) e" C0 * + C , 

for some positive constants cq,Co independent of initial data and time. The proof 
of Theorem 13.41 is now complete. 

5. Appendix 

We will consider a more general problem than (|1.7[) - (|1.8p by taking fi(s) = 
f (s) — As, g% (s) = g (s) — 7s, provided that A, 7 > are sufficiently large, and 

f(s)>-c f ,g(s)>-c g ,VseR, 

for some positive constants c/, c g . If / (s) ~ Ct \s\ p s and g (s) ~ C g \s\ q s, as \s\ — ► 
00, for p, q > 1 and some positive constants Ct,C g , it is well-known [2 3) that 
problem 

(5.1) d t u - vAu + f (u) - Ait = 0, in Q x (0, +00), 
subject to the dynamic condition 

(5.2) d t u + vd n u + g (u) — ju = 0, on T x (0, 00) , 
and initial condition 

(5.3) u\ t=0 = u o in 
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possesses a finite dimensional global attractor .4dyn which is bounded in H 2 (ft) (~1 
X°°. 

Let u* be a constant (hyperbolic) equilbrium for the system (|5.1j) - ([5.2[) (see |22[ 
Section 31). We linearize (I5.1I)-(I5.2I) around it*. We obtain 



(5.4) d t u = vAu - (f (it*) - A) u, infix (0, +oo), 
subject to the dynamic condition 

(5.5) dtu = —vd n u — (^g (u*) — 7^ u, on V x (0, 00) . 

We aim to better understand the nature of the (invariant) unstable eigenspace E u 
which corresponds to the following (matrix) operator 

h(u*)w=( uA ™- f ' ( ; u ; ) ™ + Xw ),w=( w 

with a (L (u*)) C {C : C > 0}- We note that (L (u„) , dom (L (w*))) is self-adjoint on 
X 2 with spectrum contained in (—00, C\ n ] , for some C\ n > which depends only 
on /, <7, A and 7 (see, e.g., [33] and references therein). Next, let {ipj ( x )}j et $ t> e 
an orthonormal basis in X 2 consisting of eigenfunctions of the (positive) Wentzell 
Laplacian Aw operator (see [22j Theorem 5.1]) 

(5.6) A w <Pj = &j<Pj, J e No, ifj € dom (A w ) n C (H) 
such that 

= A < Ai < A 2 < ... < Aj < A J+ i < .... ->• +00. 
We shall seek for eigenvectors Wj — (™ j ) £ X 2 , of the form Wj (x) = (fj (x)pj, 
Pj £ R, satisfying equation 

(5.7) L (u*) W 3 = CjWj, W 3 e dom (L («,)) := dom (A w ) . 

Note that for Wj £ dom (L (u*)) C H 1 (fi) x L 2 (r) , the trace of Wj makes sense as 
an element of H 1 / 2 (T). Substituting such Wj into (|5.7p . taking into account (|5.6|) 
and the fact that 

(-/' («*) + A)w,- 



(-.9' (u*)+j)Wj 



L (u.) = -uA w Wj + H\^Wj, U x ,jWj := 
we obtain the equation 

' -/'(«,) + A 



(5.8) | , I ■ Ux - .p, ( jPj , 11a { 

V " — 9 \ u * ) + 7 

A nonzero Pj exists if C = Cj i s a root of the equation 

(5.9) det {-vAjl + IL\, 7 - CO = 0, C > 0. 

When ^ = 0, this equation has at least one root £ > provided that at least one of 
A and 7 is sufficiently large, i.e., either A > / (u„) or 7 > g (u*) (in fact the roots 
are £ = A — / (m,) and £ = 7 — 5 (it*), respectively). Therefore, there exists 5 > 
such that when i/Aj < (5, the equation (|5.9|) has a root Cj (L) = (^) with Cj > 0. 
Therefore, to any such root Cj, we can assign a nontrivial pj, which is a solution of 
(|5.8p . and thus an eigenvector Wj. Let us now compute how many j's satisfy the 
inequality vAj < 6. When N > 3, the asymptotic behavior of Aj is 

(5.10) Aj ~ C s (T)j 1 ^ N ~ 1 '> as j 00 
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(see [22l Theorem 5.4]). The inequality i/Aj < 8 certainly holds when 



(5.ii) i<j<c x „S n - 1 (c 8 (p)v) 1 - N = 




, for N > 3, 



where the positive constants CA. 7 , C a depend only on A, 7 and N. 

Remark 5.1. Note that the number of unstable mode solutions to (|5.4p - (|5.5l) obeys 
the same relation (I5.11[) even when / = and A = in (|5.1|) (i.e., the dynamics of 
u inside the bulk SI is strictly linear). Finally, we note that both G\ 7 , C x — > +00 
if either 7 — > +00 or A — > +00 (cf. also [32~J Section 3]). In this case the instability 
index of u* is 

N + («.) ~ C Xn \T\ (- , N > 3. 
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